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Multi-dimensional fission model with a complex absorbing potential 

Guillaume Scamps^’B and Kouichi Hagino^’^’0 

^Department of Physics, Tohoku University, Sendai 980-8578, Japan 
^Research Center for Electron Photon Science, Tohoku University, 1-2-1 Mikamine, Sendai 982-0826, Japan 

We study the dynamics of multi-dimensional quantum tunneling by introducing a complex ab¬ 
sorbing potential to a two-dimensional model for spontaneous fission. We fist diagonalize the Hamil¬ 
tonian with the complex potential to determine a resonance state as well as its life-time. We then 
solve the time-dependent Schrodinger equation with such basis in order to investigate the tunneling 
path. We compare this method with the semi-classical method for multi-dimensional tunneling with 
imaginary time. A good agreement is found both for the life-time and for the tunneling path. 
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I. INTRODUCTION 

A spontaneous fission is a typical example of multi¬ 
dimensional quantum tunneling, and its description has 
remained a challenge in nuclear theory. The first step for 
any calculation is to construct a potential energy surface 
in a multi-dimensional space. A standard approach for 
this is the macroscopic-microscopic method, with which 
the potential surface is constructed using the shell cor¬ 
rection method of Strutinsky EH- In recent years, a 
microscopic description based on a self-consistent mean- 
field theory has also been attempted [M3- Even though 
a microscopic understanding of fission phenomena is im¬ 
portant, there have still been many open problems to be 
solved. For instance, the choice of relevant degrees of 
freedom is still under discussion Q and a large uncer¬ 
tainty may arise from a choice of energy functional M- 
Moreover, a difficulty in constrained mean-field calcula¬ 
tions has also been pointed out H, [13 ■ 

In order to calculate the fission life-time, most of the 
calculations, both with the macroscopic and the micro¬ 
scopic approaches, rely on the semi-classical approxima¬ 
tion. That is, one often searches the least action path in 
a multi-dimensional space flan or equivalently solves 
the Newtonian equations with the inverted potential [13- 

[13- 

In this paper, we investigate this problem from a dif¬ 
ferent perspective. That is, we solve the time-dependent 
Schrodinger equation (TDSE) and monitor the time evo¬ 
lution of wave function in a fully quantum mechanical 
manner. This method provides a good intuitive descrip¬ 
tion for particle decays, and has been applied to systems 
where the decay width of a resonance state is about the 
same order of magnitude as the resonance energy [13-[13 ■ 

In the previous applications of this method, because 
the TDSE has been solved by a numerical integration, the 
time step had to be small with respect to the typical time 
scale of the process. The evolution was then restricted to 


the small time regime. In this paper, we propose a new 
method introducing a complex absorbing potential in the 
exterior region. The complex absorbing potential has 
been used in time-dependent calculations [2§ - l33| in or¬ 
der to absorb the wave function at the boundary so as to 
avoid the reflexion that would otherwise perturb the dy¬ 
namics. For the decay problems, the absorbing potential 
simulates the outgoing wave boundary condition, which 
is imposed when one constructs a Gamow state. Notice 
that this method is intimately related to the so called 
complex absorbing potential method (the CAP method), 
which has been developed to compute resonances states 
in atomic physics |33l - l35l| as well as in nuclear physics 
[ 334 ^ . even though we do not take the limit of vanish¬ 
ing complex absorbing potential. 

The paper is organized as follows, In Sec. m we present 
the formalism for time-dependent calculations with a 
complex absorbing potential, which describes quantum 
mechanically multi-dimensional tunneling decay prob¬ 
lems. In Sec. imi we apply this method to a simple 
two-dimensional model for spontaneous fission. We com¬ 
pare the results with those in the semi-classical approxi¬ 
mation both for the one- and two-dimensional problems. 
We then summarize the paper in Sec. |TVl 


II. FORMALISM 

In order to investigate the multi-dimensional tunneling 
problem, we solve the TDSE in a finite box. To this end, 
we add an imaginary potential to the Hamiltonian H, 
that is, H'=H -\- iW{r). The imaginary potential iW{r) 
absorbs the outgoing flux, and should be applied only at 
the edge of the box in order not to perturb the physical 
behavior of the decay process. This is effectively equiv¬ 
alent to imposing the outgoing wave boundary condition 
for resonance states. The TDSE can be integrated as, 

|vI/(t))=e-^*^>o), (1) 


where |^o) is the initial wave function at t=0. 
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formed by the left and right eigenfunctions of the Hamil¬ 
tonian H', 

H'\^:) = E,\^:) and = E^l (2) 


with the barrier height of 


Vb 


1 / 2Rb 
6 V3MH2 



muj‘^ J 


3 


( 8 ) 


Notice that the eigenvalues Ei are complex since the 
Hamiltonian E[' is non-Hermitian. The bi-orthogonal ba¬ 
sis forms the completeness relation as = 1, 

which leads to the simple evolution, 

|vl,(t))=^e-i‘^‘(^'|vI/o)|</^r). (3) 

i 

An advantage of this method is that the evolution of the 
system can be followed for a very long time. This is 
particularly suitable for a tunneling process where the 
life-time is several orders of magnitudes longer than the 
characteristic time scale of the system. 

Among the eigenstates of El', we identify the state \ip'[) 
which has the smallest value of the imaginary part of 
eigenenergy, Ei=El — iTi/2, with the physical resonance 
state. The real part of energy, E^ corresponds to the res¬ 
onance energy while the imaginary part T i corresponds to 
the resonance width. In fact, it is straightforward with 
the TDSE to show that this state has the life-time of 
Ti=h/Ti. In this method, the details of the initial wave 
function j'l'o) is unimportant as long as it has an appre¬ 
ciable overlap with the resonance wave function, \^pl). 


When the intrinsic degree of freedom ^ is neglected, 
the saddle is at Rs=2Rb/3 with the height of 
Vb=2Mn^Rl/27 [12. 

In the calculations presented below, we use the same 
parameters as those in Refs. [l^l, 112 except for 
Rb, which we vary to study the tunneling in the po¬ 
tential with different barrier heights. Those param¬ 
eters were determined in order to mimic the sym¬ 
metric fission of with a coupling to the beta 

vibration. The parameters are then taken to be 
Hn=hu}=0.97 MeV, g^=Mmn^M=23 AMn/^ 
and m=3AMNRQ/{8Tr), where Mjy is the nucleon mass, 
A=234 is the atomic number of the nucleus, and 
i?o=l-2A^/^ fm is the equivalent sharp radius (notice that 
the vertical axis in Fig. 4 in Ref. is actually Rq^, 
rather than ^ itself). In the actual calculations, we mod¬ 
ify the total potential to a constant value in the outer 
region in order to avoid the divergence of the potential, 
that is, V{R,^) max(I4(i?,f), —3). 

B. 1-Dimensional Problem 


III. RESULTS 
A. Model Hamiltonian 

We now apply the formalism presented in the previous 
section to a fission problem and compare the results with 
those in the semi-classical approximation. To this end, 
we em ploy a two-dimensional fission model considered in 
Refs. [ 42 , 13 . This model consists of the elongation R 
between the two fission fragments and an intrinsic degree 
of freedom ^ coupled to the elongation. The Hamiltonian 
then reads. 


Before we discuss the tunneling dynamics in the two- 
dimensional space, let us first solve the problem in one 
dimension neglecting the ^ degree of freedom. In the 
one-dimensional problem, the tunneling path is trivial, 
and the decay life-time is obtained in the semi-classical 
approximation as, 

(9) 


S= y/2M{U{R) - Eq) dR, (10) 

jRo 


9 ^ 1 2^2 DA 


The potential U(R) is chosen to be 

U{R) = fl- ^ 

2 V Rb 


(4) 


(5) 


in order to form a barrier. The total potential, 
V{R,^)=U{R) + E has a saddle at 


Ra 

is 


2Rb 

3MER 



muj'^ ) ’ 


2gRb 

3MmH2a;2 



mu}^ J ’ 


( 6 ) 

(7) 


where Eq is the energy of the decay state, and Rq and 
Ri are the inner and the outer turning points, respec¬ 
tively, satisfying U{Rq)=U{Ri)=Eq. For a cubic poten¬ 
tial given by Fq. (ISl). the semi-classical formula can also 
be transformed to [20, Ell , 


with 


T = 


1 

o 




^2MU{R) dR, 


( 11 ) 


( 12 ) 


where [/(R'l) = C/(0) = 0. 

In order to solve the same problem quantum mechan¬ 
ically, we first determine the initial wave function Id/o) 
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by modifying the potential U{R) so that the modihed 
potential has a bound state. To this end, we replace 
the potential outside the barrier with a constant value 
1^ , which we take 1 MeV as shown by the dashed 
line in Fig. [T] (a). The corresponding initial wave func¬ 
tion is plotted by the dashed line in the lower panel of 

Fig.m 
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FIG. 1: (Color online) (a) The one dimensional potential 
U{R) with = 6 fm (the red solid line). The figure also 
shows the modified potential for the initial wave function 
(the blue dashed line) and the absorbing potential W (R) (the 
green dotted line), (b) The resultant initial wave function 
(the blue dashed line) and the resonance wave function (the 
red solid line). 

For the complex absorbing potential, we employ the 
shifted polynomial function , 


W{R) = Wo{R - Rafe{R - Ra), 


(13) 


with Wo=~0.387 MeV and i?a=14 fm. We have con¬ 
firmed that the results do not significantly change even 
if we vary the value of Wq and Ra as well as the size of 
the box. 

With this complex potential, we then determine the 
ensemble of eigenstates and eigenvalues of the Hamilto¬ 
nian H' defined by Eq. ([2]) . It should be mentioned here 
that due to the large difference in the order of magnitude 
between the real and imaginary parts of the resonance 
energy, it is numerical necessary to use the quadrupole 
precision in the program. With this prescription physi¬ 
cal quantities can be calculated up to about 34 decimal 
digits. This allows us to calculate a life-time of the or¬ 
der of billion of years when the characteristic time of the 
system is of the order of 10“^^ s. Another important pa¬ 
rameter in the calculations is the lattice mesh size Ai?, 
that has to be small enough in order to describe correctly 


the tunneling wave function. In our calculations, we take 
Ai?=0.25 fm with the finite difference formula with 9 
points for the second derivative in the kinetic energy op¬ 
erator. 

To select the physical resonance wave function among 
the eigenstates of H', we take the lowest energy state 
\ip\) which has the maximum overlap with the initial 
state. The resulting wave function is shown in Fig. [1] 
(b) by the solid line. We see that the wave function is 
smoothly damped up to a factor of 10^° before reaching 
the reflecting edge of the lattice at i?=54 fm. The decay 
width can then be read off from the imaginary part of 
the eigenenergy. 





R [fm] 

FIG. 2: (Color online) The time-evolution of the square of the 
wave function during the decay process. The time is indicated 
in unit of fm/c. 

In order to have a more intuitive picture on the de¬ 
cay process, we compute the time-evolution of the wave 
function according to Eq. The probability density at 
different times obtained with a larger value of Ra (i?a=54 
fm) is shown in Fig. [2j At t=0, the confining in the mod¬ 
ified potential is suddenly removed, and the initial wave 
function is coupled to the continuum. In the first instant 
of the dynamics, one can see the emitted wave going out¬ 
ward the potential barrier. After about t=10^ fm/c, the 
emission becomes stationary, and the wave packet has the 
same shape as the metastable wave function constructed 
by diagonalizing the Hamiltonian (see the solid line in 
Fig. [11(b)). At later time, the wave packet is absorbed 
exponentially keeping the same spatial shape. 

The decay width can be calculated with the TDSE by 
computing the survival probability defined as 


V{t) = 


rAi. 


\^{R,t)\'^dR, 


(14) 


as a function of time, where Aum is taken outside the bar¬ 
rier. The survival probability obtained with i?iim=20 fm 
is plotted in Fig. [3] for the potential with i?{,=ll fm. 
We find that both the TDSE method and the imagi¬ 
nary part of the eigenenergy of H' yield 1/t=5.76 x 10“® 
(1/year) while the semi-classical approximation yields 
1/t=5.69 X 10“® (1/year). The agreement between the 
quantal and the semi-classical calculations is rather good 
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for this parameter set, partly because the barrier is high 
and the multiple reflections under the barrier can be ne¬ 
glected. 
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FIG. 3: (Color online) The survival probability as a function 
of time computed dynamically using the formula m with 
I?iim=20 fm (the blue crosses). The Rb parameter in the po¬ 
tential is taken to be Rb=ll fm. The survival probability is 
plotted both in the linear scale (the upper panel) and in the 
logarithmic scale (the lower panel). For a comparison, the 
figure also shows the survival probabilities for an exponential 
decay with the decay width obtained from the semi-classical 
approximation given by Eq. ([9]) (the red solid line). 


C. 2-Dimensional Problem 

We now discuss the tunneling dynamics in the two- 
dimensional surface. To this end, we diagonalize the 
Hamiltonian H' in 2 dimensions in a lattice of dimen¬ 
sion R € [—4.75 fm : 16 fm] and ^ G [—7.25 : 5.5] with a 
mesh size of Ai?=0.25 fm and A^=0.25. The imaginary 
potential is implemented with the same expression as Eq. 
(fT^ with Wo=—0.1 MeV and i?a=ll fm- The reference 
wave function to be used to select the physical resonance 
state is obtained in a similar manner as in the previous 
subsection. The wave function for the metastable state 
is plotted in Fig. |4| for the choice of i?t,=6 fm, together 
with the total potential V{R, ^) in the contour lines. One 
can see that this wave function is well confined inside the 
barrier with a small component outside due to the quan¬ 
tum tunneling, as in the one-dimensional case shown in 
Fig. [H 

Notice that this wave function corresponds to the wave 
function after the decay becomes stationary. The pre- 
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FIG. 4: (Golor online) The wave function for the meta-stable 
state plotted in the logarithmic scale. This state corresponds 
to an eigenstate of the Hamiltonian H' for the potential with 
Rb~G fm and with the imaginary potential given by Eq. ca) 
with Wo=—0.1 MeV and 7?a=ll fm. The total potential 
V {R, is also plotted by the contour lines. 

stationary decay at the first instant of the decay process 
is shown in Fig. [5] following the time evolution from the 
reference wave function as the initial state. One can see 
that the flow occurs in a small region in the potential 
energy surface. 
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FIG. 5: (Color online) The time evolution of the two- 
dimensional wave function at time t=0, t=250 fm/c, t=500 
fm/c and t=1250 fm/c. For a comparison, the semi-classical 
tunneling path is also shown with the dashed line. 

We compare the quantum mechanical flow with the 
semi-classical tunneling path. To this end, we follow 
the method described in Refs. [20l - [^ . That is, the 
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semi-classical tunneling path is determined by solving 
the classical equations of motion in an inverted poten¬ 
tial, V{R,^) —V{R,^). For the initial condition, 

one can use R{0)=e^ 2 MU cos0, sin0, and 

i?(0) = ^(0)=0, where e is a small number. By searching 
with 9 from 0 to 27r, one finds a special angle 9 for which 
the classical path reaches the equi-potential surface (and 
bounces back to the origin if one continues to follow the 
time evolution of the path), while for other values of 9 
the classical path is reflected before it reaches the equi- 
potential surface [20l - l^ li^ . . This special path is re¬ 
ferred to as the escape path, and plays an important role 
in the semi-classical theory of multi-dimensional quan¬ 
tum tunneling. The escape path so obtained is denoted 
by the dashed line in Fig. [5] (there is only one escape path 
for the potential considered in this paper). It is remark¬ 
able that the quantum mechanical time evolution almost 
follows the semi-classical path. We have confirmed that 
the tunneling path obtained by minimizing the classical 
action with the algorithm in Ref. [d^ provides the same 
result. 

To make a further comparison of the tunneling path, 
we compute the flux, j={jR, j{), from the resonance wave 
function as. 


and 


jR 


ft ( r* M 

2iM V ^ dR 










(15) 


(16) 


The flux is shown in Fig. Inland is compared to the semi- 
classical escape path. The quantum flow is a collection 
of all the trajectories connecting the region around the 
origin with the continuum region. We see that the main 
quantum mechanical flow is systematically parallel to the 
semi-classic trajectory. The semi-classical path can thus 
be regarded as the mean trajectory of the quantum flow, 
where as the quantum flow takes into account the quan¬ 
tum mechanical fluctuation around the classical trajec¬ 
tory. 

In addition to the agreement for the tunneling path, 
we also compare in Fig. [7] the resulting fission life-time. 
Varying the Rb parameter in the potential U (i?), we ob¬ 
tain a range of fission life-time from the characteristic 
time of the system to a life-time longer than the actual 
life-time of the nudeug xhe quantum mechanical 
life-time is computed from the imaginary part of the res¬ 
onance energy, while the life-time in the semi-classical 
approximation is evaluated using the formula Eq. m 
with the action integral evaluated along the escape path 
V in the imaginary time, 


S'= J 


dr. 


(17) 


The latter is equivalent to approximating the pre¬ 
exponential factor by that for a one-dimensional problem 



R [fm] 

FIG. 6: (Color online) A comparison of the quantum mechan¬ 
ical flux computed with Eqs. (1151) and (1161) (the red arrows) 
with the semi-classical escape path (the blue dashed line). 
The two-dimensional potential energy surface is also shown 
by the black contour lines. 


with a cubic potential along the escape path [d^. In the 
range shown in the hgure, a good agreement is found be¬ 
tween the two methods, the maximum deviation being 
up to about 20%. Evidently, the semi-classical method 
provides a good approximation to the multi-dimensional 
tunneling problem for this Hamiltonian. 



Rb [fm] 


FIG. 7: (Color online) The fission life-time as a function of the 
Rb parameter in the potential U{R). The life-time computed 
quantum mechanically is denoted by the blue crosses, while 
that evaluated with the semi-classical method using Eq. 113 
is denoted by the red line. 


IV. SUMMARY 

We have presented a full-quantum method to study a 
decay life-time as well as the tunneling dynamics. To this 
end, we have introduced a complex absorbing potential, 
which is effectively equivalent to the outgoing boundary 
condition, and constructed the bi-orthogonal basis. We 
have shown that this method provides a good tool to fol¬ 
low the time evolution of a system over a very long time. 
This enables one to compute the decay life-time from a 
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few fm/c to the order of billion of years. A comparison 
with the semi-classical approximation has shown a good 
agreement between the two methods. It has been shown 
that the average of fissioning flux in a multi-dimensional 
plane corresponds to the semi-classical tunneling path. 
For the decay life-time, the two methods yield similar 
values to one another, where the maximum difference is 
only about by 20%. 

The quantum mechanical method discussed in this 
paper provides a good alternative to the semi-classical 
method. It provides an intuitive picture of multi¬ 
dimensional quantum tunneling, including the quantum 
mechanical fluctuation of the classical path. It will also 
provide a convenient method when a bifurcation of tun¬ 
neling path is important, e.g., in the presence of a com¬ 
petition of several fission modes. We plan to apply 


this method, both for dynamics and a determination 
of resonance state, to more realistic systems using the 
constrained mean field theory for the potential energy 
surface and/or the time-dependent generator coordinate 
method We will report it in a separate pub¬ 

lication. 
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